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ABSTRACT 

Observations of quasar absorption line systems reveal that the z = 3 intergalactic 
medium (IGM) is polluted by heavy elements down to HI optical depths tri <C 10. 
What is not yet clear, however, is what fraction of the volume needs to be enriched 
by metals and whether it suffices to enrich only regions close to galaxies in order to 
reproduce the observations. We use gas density fields derived from large cosmological 
simulations, together with synthetic quasar spectra and imposed, model metal distri- 
butions to investigate what enrichment patterns can reproduce the observed median 
optical depth of C IV as a function of tri- Our models can only satisfy the observa- 
tional constraints if the z = 3 IGM was primarily enriched by galaxies that reside in 
low-mass (m to t < 1O 1O M0) haloes that can eject metals out to distances > 10 2 kpc. 
Galaxies in more massive haloes cannot possibly account for the observations as they 
are too rare for their outflows to cover a sufficiently large fraction of the volume. 
Galaxies need to enrich gas out to distances that are much greater than the virial 
radii of their host haloes. Assuming the metals to be well mixed on small scales, 
our modeling requires that the fractions of the simulated volume and baryonic mass 
that are polluted with metals are, respectively, > 10% and > 50% in order to match 
observations. 

Key words: galaxies: formation — intergalactic medium — quasars: absorption lines 
— methods: N-body simulations 



1 INTRODUCTION 

Studies of quasar absorption spectra indicate that at 
redshift z > 2 most of the bar yons in the Uni- 
verse reside outside o f galax i es (e .g. iRauch et al ] Il997l : 
IWeinberg et al.1 Il997l : ISchavd l200lh in the gas that is 
observable through HI Lya absorption and that pollu- 
' tion of the low-density intergalactic medium (IGM) by 
heavy elem e nts i s wide- spread (e.g . iTvtler et al. [l995: 
Cowie et all Il995l: ICowie fe Songailal Il998l: lEllison et all 
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Aracil et al.l 120041: ISimcoe et al. 2004; IPieri fe Haehneltl 
20041 : ISongaildl2005l : lAg"uirre et al.1 120081 : IPieri et al.ll2010h . 
The very fact that heavy elements are able to make their 
way from galaxies out into intergalactic space, has far 
reaching implications for the galaxy formation process. We 
know that metals are only produced in the high-density 
environments where star- format ion occurs, and in order 
to reach the low-density IGM, they must have been re- 
moved from these regions. Energetic feedback processes that 
drive gas out of galaxies are believed to be the most im- 
portant way in which the IGM is polluted by metals, as 
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suggested by both observational (e.g . IPettini et al.l l200ll : 
S ha p lev et alj 120031; ISteidel et al.ll2010l) and theoretical (e.g . 
l Aguirre et al.1 l2001bl: jTheuns et al.1 120021: ICen et al.1 120051 : 
iQppenheimer fc D ave 2006: iTornatore et alj feoiQ) studies. 

In spite of the observational evidence that metals have 
polluted the high-redshift IGM, several issues remain un- 
clear. What fraction of the volume of the IGM, and hence of 
the Universe, is polluted by metals? Out to what distance 
do galaxies need to enrich the gas around them in order to 
reproduce observations? Can observed galaxies do the job 
or are fainter galaxies the main culprit? 

The metal distribution in the high-redshift IGM 
has been investigated by many authors using self- 
consi ste nt, hydrodyn am ical si mulations (|Theuns et akl 
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IQppenheimer &; Davef 

Qppenheimer &; Pavel 12008: _ 

120101 : IWiersma et al.1 l2009bL l20ld : ICen fc Chisaril l201lh . 
The interpretation of these simulation results is, however, 
complicated by their complex nature and by the crudeness 
of, and freedom provided by, the required subgrid models. 
In addition, their computational expense prohibits compre- 
hensive explorations of parameter space. Observationally, 
IPieri fc Haehneltl (2004) use four quasar spectra along with 
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a suite of synthetic quasar spectra and find that the lower 
limit for the volume filling factor of O VI is > 4%. 

Simple models therefore represent a useful comple- 
ment to full-blown simulation studies. To this end various 
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metals placed around haloes to investigate the metal dis- 
tribution in the universe. Models in which the gas is en- 
riched down to varying densities have been use d to constrain 
the volume filling factor of e nriched gas (e.g. ISchave et al.l 
2003; Pieri & Haehnelt 2004). However, as discussed by e.g. 



Schave &; Aguirrd (J2005|), the inferred filling factors could 
be misleading if the metals are poorl y mixed, as oberva tions 
suggest to be the case on both la rge (Sch ave et al.| [2003) and 
small scales ([Schave et al.ll2007h . 

In the present work we combine toy models for the 
metal distribution with a large, cosmological hydrodynami- 
cal simulation to investigate what halo masses could host the 
sources of the observed intergalactic metals, out to what dis- 
tances the galaxies need to enrich the gas, and what (large- 
scale) volume filling factor of heavy elements is required in 
order to reproduce the observed metal distribution as probed 
through absorption lines in the spectra of quasars. In par- 
ticular, we compare model predictions to the observed z — 3 
relation between the median optical depth of C I V as a func- 
tion of the H I optical depth ([Schave et al.ll2003h . We choose 
to restrict our analysis to the median optical depth because 
it allows us to compare with published results and because 
it provides a simple measure of how far away from galaxies 
the IGM is being enriched. 

We achieve this by extracting synthetic absorption spec- 
tra after imposing simple metal distributions in which all 
haloes above a given mass cut are allowed to enrich the 
IGM spherically out to a fixed radius. This allows us to link 
a metal distribution with a well-defined mass and volume 
filling factor, via cosmological gas density and temperature 
distributions, to the observations and so to determine which 
of the possible metal distributions are capable of reproduc- 
ing the observations. We will show that the IGM must have 
been primarily enriched by galaxies that reside in low-mass 
(ratot < 10 10 M©) haloes and are capable of driving gas 
out to distances > 10 2 kpc. Assuming the metals to be well 
mixed on small scales, the fractions of the volume and bary- 
onic mass that are polluted with metals are, respectively, 
> 10% and > 50% in all models that are capable of match- 
ing the observations. 

This paper is organised as follows. We first introduce 
our simulation set (Sec. 12. ip and the methods we use to 
distribute the metals (Sec. 12.2]) , then, in Sec.[3l we describe 
our results and in Sec. [H we summarise our findings and 
conclude. 



2 METHOD 
2.1 Simulations 

The simulation analysed in this study is one of the cos- 
mological, hydrodynamical simulations that comprise the 



Table 1. Simulation parameters. From left to right: Simulation 
identifier; comoving box size (I/box); number of particles of both 
dark matter and gas; mass of dark matter particles (ttidm)? and 
initial mass of gas particles (m ga s)- 



Simulation Lb ox 
(Mpc//i) 



N 



WDM 

(M©) 



(M©) 



L025N512 
L025N256 
L012N512 
L012N256 



25.0 
25.0 
12.5 
12.5 



512 3 
256 3 
512 3 
256 3 



8.68 x 10 6 
6.95 x 10 7 
1.09 x 10 6 
8.68 x 10 6 



1.85 x 10 6 
1.48 x 10 7 
2.31 x 10 5 
1.85 x 10 6 



Overwhelmingly Large Simulations (OWLS ) project and 
is described in detail in Schav e et al ] (|2Q10h . Briefly, the 
simulation was run using a significantly extended version 
of the parallel PMTree-Smoothed Particle Hydr odynam- 
ics ( SPH) code gadget in (last described in ISpringell 
2005), a Lagrangian code used to calculate gravitational 
and hydrodynamic forces on a particle by particle ba- 
sis. The simulations track star formation, supernova feed- 
bac k, radiative cooling an d chemodynamics, as described 
i n ISchave fe Dalla Vecchial (I2008T ). i Dalla Vecchia & Schavel 
(|2008MWiersma et al.1 (|2009ah and lWiersma et al.1 ()2Q09bh 7 
respectiv ely. This physical m odel is denoted as the REF 
model in Sc have et al.l (|2010h , and is used in all of the sim- 
ulations analysed in this paper. 

For the purposes of this work, our prescriptions for ra- 
diative cooling and reionization are the most important as- 
pects of the model, as the thermal state of the IGM depends 
on them. In brief, we c alculate radiati v e cooli ng and heat- 
ing using the tables of Wiersma et al ] (|2009ah . which con- 
tain net cooling rates (calculated u sing the code cloudy, 
last described in iFerland et ai1ll998h as a function of den- 
sity, temperature and redshift for each of the 11 elements 
hydrogen, helium, carbon, nitrogen, oxygen, neon, magne- 
sium, silicon, sulphur, calcium and iron, computed under the 
assu mption of ionization eq uilibrium and in the presence of 
the lHaardt &; Madau] ([200 lh model for the uniform, evolv- 
ing meta-galactic UV and X-ray radiation field from galax- 
ies and quasars as well as the cosmic microwave background. 
The si mulations model hydroge n reionization by switching 
on the lHaardt & Madaul (|2Q0lh background at z = 9. He- 
lium reionization is modelled by heating the gas by a to- 
tal amount of 2 eV per atom. This heating takes place at 
z = 3.5, with the heating spread in redshift with a Gaus- 
sian filter with cr(z) = 0.5. This reionization prescription 
used in these simulations matches the t emperature h i story o f 
the IGM infe rred from observations by ISchave et al.l (|2000bh 
(see Fig. 1 of I Wiersma et alJl2009bh . 

The simulations used in this paper are summarised 
in Table [T] All assume a flat ACDM cosmology with 
the cosmological parameters: {n m , ^b, ^a, 08, n s , h} = 
{0.238, 0.0418, 0.762, 0.74, .951, 0.73), as deter mined from 
the WMAP 3-year data dSpergel et al.1 l2007h - The sim- 
ulations from which we derive the bulk of our results 
are L025N512 and L012N512. Two additional simulations, 
L025N256 and L012N256, are used to independently assess 
the effects of simulation box size and numerical resolution. 

Although each of the simulations was run to z — 2 
we restrict our analysis to the z = 3 simulation snapshots, 
approximately corresponding to the median redshift of the 
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observational sample that we compare to. Our analysis de- 
pends on the identification of the masses and locations of 
gravitationally bound dark matter haloes, which are identi- 
fied using the spherical o verdensity criterion i mplemented in 
the SubFind algorithm (|S pringel et al.ll200ll ). Halo proper- 
ties quoted in this paper are defined with respect to spheres 
with radius r2oo and mass 772200, centred on the potential 
minimum of each identified halo, defined so that they con- 
tain a mean internal density equal to 200 times the critical 
density of the Universe at the redshift we are considering. 

We note that although we will present results only for 
the reference implementation of the subgrid physics mod- 
ules, we have repeated the analysis for a range of physics 
implementations. This is important because the different 
physics prescriptions can affect the density, temperature 
and velocity fields of the absorbing gas, changing the pre- 
dicted ion abundances and optical depths. We find, however, 
that using either a simulation with strong, AGN feedback 
(AGN.L025N512 in the OWLS nomenclature), or a simu- 
lation that neglects both supernova feedback and cooling 
through metal lines (NOSN_NOZCOOL_L025N512 in the 
OWLS nomenclature) has a negligible effect on our results 
or conclusions. 



2.2 Imposed metal distributions 

In order to predict a synthetic rciv — thi relation, we re- 
quire knowledge of both the distribution of metals and the 
physical state of the absorbing gas. In the simulation, gas 
metallicities are tracked self-consistently, but in the present 
study we do not make use of this information. We instead 
assume that all haloes with a total mass greater than m z are 
able to enrich the surrounding gas out to a proper distance 
r z to a metallicity Z, and that outside of these spheres the 
metallicity of the IGM is zero. Our model for the intergalac- 
tic metal distribution is therefore completely specified by 
three parameters: ra z , r z and Z. As we will see, the parame- 
ters r z and m z determine the shape of the relation between 
rciv and thi, which is the primary focus of this paper. The 
metallicity changes only the normalisation, with rciv oc Z, 
so we simply scale Z in each run to match the normalisation 
of the observed rciv — thi relation at log 10 (rHi) = 2.5, the 
largest optical depth probed by the observations. 

At this point we must note one caveat: we have imposed 
metal distributions on to already completed simulations, the 
models are not fully self- consistent in that they do not in- 
clude the effect of the winds that carry the metals on the 
density and temperature structure of the gas and in that 
they do not include the effect of the metals on the cooling 
rates. The last point is, however, not a major concern as we 
will show that the metallicities required to match the obser- 
vations are sufficiently low (typically 10 -3 — 10 -2 Z©; Fig. [2]), 
that metals do not significa ntly change the gas cooling rates 
(e.g. IWiersma et"aIll2009eJ ). 

While the hydrodynamical simulations that underlie 
our models did include winds, these winds fall short of be- 
ing able to account for the observed C IV at low thi, as we 
will show elsewhere. This failure is actually consistent with 
our results. Although the simulations have sufficient reso- 
lution to identify dark matter haloes to very low masses 
(corresponding to ~ 10 2 dark matter particles), their fi- 
nite resolution does cause us to strongly underestimate the 



star formation rates in most of the low-mass haloes that 
we can identify. Hence, the simulations underestimate the 
number and strength of the outflows originating from the 
low-mass haloes that we claim to be responsible for the en- 
richment of the IGM. Because not all of the gas that is en- 
riched in our models was touched by winds in the underlying 
hydro simulation, we cannot exclude the possibility that a 
self-consistent simulation giving rise to a similar distribu- 
tion of bubbles would predict the enriched gas to be too 
hot to be visible in C IV. Reassuringly, we find, as noted 
above, that post-processing simulations without winds or 
with much stronger winds leads to identical conclusions. 

We note that we expect the heating effect of outflows 
from low-mass galaxies to be smaller than those from the 
more massive galaxies that our simulations do include. This 
is because there are already many low-mass galaxies at high 
redshift, giving the gas more time to cool, and because they 
are observed to drive winds of velocit ies < 10 2 km/s (e.g. 
ISchwartz &; Martini [20041 : lMartinll2005l ) (this velocity corre- 
sponds to post-shock temperatures of ~ 2 x 10 5 K, assum- 
ing the gas is fully ionized and of primordial composition), 
which leaves the post-shock gas at temperatures for which 
the cooling time is muc h shorter than the age of the Universe 
(jWiersma et alJl2Q09aJ ) . 

The conclusions based on our simple models will, how- 
ever, ultimately need to be confirmed by self- consistent, hy- 
drodynamical simulations. Unfortunately, at the moment 
such simulations rely on uncertain subgrid models fo r the 
generation of winds fe.g. lDalla Vecchia &; Sch ave 2008) and 
they lack the resolution required to model outflows from 
low-mass galaxies and to simu late the small-sc ale mixing 
relevant for the observations (jSchave et al.ll200/t ). 

The bulk of our results are derived from a grid of mod- 
els in which m z is varied in steps of 0.5 dex from the lowest 
mass haloes that can be robustly identified in the highest 
resolution simulation, m z — 1O 8 M0, which are the least 
massive haloes that ar e expected to be able to produce 
stars after reionization jEfstathioulll992l : lOuinn et all fl 996: 
iThoul fc Weinberdll996h . up to m* = 10 11 M . 

Note that a mass of 10 11 Mq is small compared with 
the total masses inferred fo r obser ved galaxies at z = 3. For 
example, lAdelberger et al.l (|2005l ) find that Lyman-break 
galaxies reside in haloes of mass ~ 10 12 M o . However, as 
we will show, such high-mass galaxies are unimportant for 
the enrichment of the IGM. 

The parameter r z is changed in factors of 2 from 
31.25 kpc to 500 kpc. In addition, we investigate a set of 
runs in which haloes in the fiducial simulation are allowed 
to enrich the IGM out to a fixed multiple of their virial ra- 
dius, r2oo- Galactic winds with velocities up to 400-600 km/s 
are frequently detected in starburst galaxies through the 
gas absorption line s that are blue-shif t ed relative to their 
host galaxies (e.g . IPettini et all l200ll : ISteidel et all l20ld : 
iRakic et af1l201lh . If we assume that winds were ejected 
from galaxies at high redshift (z ^> 3) and that their veloc- 
ities do not decrease with time, then by z — 3 galaxies can 
enrich out to a maximum radius of 0.9 — 1.4 proper Mpc. The 
assumption of a constant, high outflow velocity and launch 
at z ^> 3 make this estimate far too optimistic, but all of 
the models that match the observations require r z to be no 
large r than 500 kp c ( which i s likely still too optimistic; see 
e.g. lAguirre et alJ e.g. l2001bl ). and are thus compatible with 
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constraints placed on the metal distribution by travel-time 
arguments. 

Optical depth distributions are calculated by firing 10 3 
randomly chosen lines of sight through the simulation vol- 
ume and calculating absorption spectra for both HI and 
C IV following the pro cedure outlined in e.g. Appendix A4 
of iTheuns et all {l998). The mean HI optical depth in our 
simulations at z — 3 is T e g — 0.388, which is consistent with 
observations (e.g. ISchave et ai1l2003l : iFaucher -Giguere et al.l 
2008). In order to match observations with HIRES on the 
Keck telescope, we convolve our spectra with a Gaussian 
line-spread function with a full-width-at-half-maximum of 
6.6 km/s, and resample our spectra onto 1 km/s pixels. We 
do not add noise to our spectra, but we have verified that 
the addition of Gaussian noise with a signal-to-noise ratio 
of greater than 25 does not significantly affect any of our 
results. We generate absorption spectra for two transitions: 
HI (1215.67A) and the CIV doublet (1548.20A, 1550.78A). 
In order to compare to the observed rciv — thi distribu- 
tion, we then bin pixels in thi and calculate the median 
tc iv corresponding to the redshifts of the pixels in each H I 
bin. 

For each run, in addition to measuring the rciv — thi 
relation, we calculate the fraction of the total mass (/ m ) and 
volume (/v) that has been enriched from 

/m T^m £~ft! W 

Here, rra and hi are the SPH particle mass and smoothing 
kernel, respectively, and the sums in the numerator of each 
fraction extend only over particles with non-zero metallicity. 
We have verified that using rrn / pi instead of /if gives nearly 
identical volume filling fractions, as expected. For the solar 
abundanc^B we use the metal mass fraction Zq = 0.0127. 



3 RESULTS 

Fig. [1] shows a thin (1 comoving Mpc/h thick) slice of the 
density field at z = 3 through the centre of the L025N512 
simulation. Each panel shows a different combination of r z 
and m z . Gas with zero metallicity is shown in grey-scale, 
while metal-enriched gas is coloured red. In addition, the 
values for fy and f m are indicated for each model. Pan- 
els showing metal distributions that are consistent with the 
observations (see below) are outlined in green. 

The fractions of the mass and volume that are enriched 
do not track each other in a simple way. The ratio fm/ fv 
is always ^ 1 because metals are only placed around col- 
lapsed structures and thus preferentially in overdense re- 
gions. For m z = 10 10 M©, changing r z from 31.25 kpc to 
500 kpc changes the ratio fm/ fv from 22.7 to 1.3, as a larger 
value of r z at fixed m z allows metals to disperse further out 
of the high-density peaks. Similarly, changing m z at fixed 
r z = 250 kpc, we find that the ratio fm/ fv rises from 2.8 for 
m z — 10 9 to 7.5 for m z — 10 11 M©, because more mas- 
sive haloes are preferentially located in higher density envi- 
ronments. Given that we understand how changing the pat- 

1 This corresponds to the value obtained using the default 
abundance set of CLOUDY (version 07.02; last described by 
iFerland et 311119981 ). 



tern of metal enrichment alters the relative mass and volume 
filling factors, we now ask how this impacts the rciv — thi 
relation. 

The curves in Fig.[2]show the relation between rciv and 
thi in the synthetic absorption spectra. Each panel corre- 
sponds to a different m z . The solid lines in panels 6, c and d 
(log 10 (ra z /M©) = 9, 10, 11, respectively) show the predicted 
rciv — thi relation from the L025N512 simulation with 
various imposed metallicity distributions. The dot-dashed 
lines in panels a and b (log 10 (ra z /M©) = 8,9, respectively) 
show the predicted tciv — thi relations for the small vol- 
ume, high-resolution simulation, L012N512. In every panel 
we compare our simulated predictions to the o bserved opti- 
cal depth pixel statistics of lSchave et al.l (l200d"h for the red- 
shift r ange 2.479 ^ z ^ 4.033, as published in I Aguirre et all 
(2005), shown as the yellow points with la error bars. The 
data come from six quasar spectra, Q0420-388, Q1425+604, 
Q2126-158, Q1422+230, Q0055-269, and Q1055+461, that 
were taken with either the Keck/HIRES or the VLT/UVES. 
A full description of the sample is given in ISchave et al.l 
(2003). In each panel Z was chosen such that the r z = 
500 kpc curve exactly matches the observations at the high- 
est value of th i and the metallicity required for this normal- 
isation is given in each panel. 

In the present work we are not aiming to reproduce the 
shape of the tciv-thi relation in detail and, indeed, would 
not necessarily expect our simple models to be capable of 
doing this. Rather, we require that the models predict a me- 
dian rciv — thi relation that is consistent with, or larger, 
than observed and that the tciv-thi relation is not steeper 
than observed. We effectively combine the two conditions by 
scaling Z in each run to match the normalisation of the ob- 
served rciv — thi relation at log 10 THi = 2.5, the largest 
optical depth probed by the observations, and then requir- 
ing the model to predict median rciv at low thi that are 
consistent with or greater than observed. 

We stress that overprediction of tciv is not a prob- 
lem at small thi because our simple models make the as- 
sumption that the metallicity inside each enriched bubble is 
constant with radius. The overestimate at low tciv could 
therefore be solved by imposing a metallicity that decreases 
with radius. On the other hand, we will assume that under- 
prediction of tc iv at low th i signals the failure of the model 
because a metallicity that increases with radius is likely un- 
physical. Such an unphysical metallicity gradient would also 
have been required if we had chosen to scale the metallicity 
to match the low thi points, because the unsuccessful, con- 
stant metallicity models predict much steeper tciv — thi 
relations than observed (see Fig. [2]). In fact, such an ap- 
proach would not be possible for most of the models that 
we rule out, because they typically predict the median tciv 
to be zero at low thi. 

In principle, we could produce significantly better fits to 
the tc iv — th i relation by considering more complex mod- 
els (e.g. including metallicity gradients, anisotropic outflows 
or scatter in metallicities) , but that is not the aim of the 
present study. Given that such models would still not be 
self-consistent if the enrichment is done in post-processing, 
it is not obvious that the use of more complex models would 
be justified. 

Before proceeding, we evaluate the effect of our simu- 
lation's finite box size and numerical resolution on our re- 
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Figure 1. Thin (1 comoving Mpc/h thick) slices of the z = 3 gas density field through the centre of the L025N512 simulation for 
models with different r z (the proper radius out to which galaxies enrich the IGM) and m z (the minimum halo mass that is responsible 
for enriching the IGM). Each density field is shown in grey-scale in regions where, for a given combination of r z and ra z , the gas is 
unenriched, and in red where metals are present. Each row of plots corresponds to a different value of m z and each column shows a 
different value of r z . The numbers at the top of each panel show what fraction of the total simulation volume (/v) and gas mass (/m) 
are enriched with metals in each model. Panels outlined in green show models that satisfy the constraints provided by the observed 
rCiv — tHi relation. All models that reproduce the observations have volume filling factors > 10%. 



suits. Firstly, the dotted curves in panels 6, c and d show 
the effect of decreasing the box size by a factor of two in 
each dimension, while keeping the numerical resolution fixed 
(L012N256 vs. L025N512). It is clear that for m z ^ 10 10 M 
our results are converged with respect to box size, but that 
for m z = 1O 12 M the 12.5Mpc//i and 25Mpc//i volumes 
predict significantly different rciv — thi relations, indicat- 
ing that for haloes of this mass the simulation results from 
the 25Mpc//i volume are not necessarily converged. Note, 
however, that allowing only haloes with masses ^ 10 n M Q 
to enrich the IGM yields rciv — thi relations that are far 
steeper than those allowed by the observations. We thus con- 
clude that our simulation boxes are sufficiently large for our 
purposes. Secondly, the long-dashed curves in panels c and 
d show the effect of degrading the simulation mass resolu- 
tion by a factor of eight while keeping the box size constant 
(L025N256 vs. L025N512). Decreasing the resolution does 



not significantly change any of the conclusions derived from 
this analysis. 

Considering first the case where only haloes with mass 
> 10 11 Mq enrich the IGM, we find that even if such objects 
were able to enrich out to a radius of 500 kpc, they would 
fall far short of being able to account for the observed rciv 
at thi < 10 because there are too few of them to enrich 
enough of the low-density gas. If we now consider a lower 
value of m z = 1O 9 M , the lower clustering strength and 
higher number density of these less massive haloes allows 
them to pollute lower density gas, even if they enrich out to 
much smaller distances. All values of r z ^ 125 kpc provide 
a good match to the observations (note that measurements 
of rciv are upper limits for log 10 (rHi) < 0.3). The inter- 
mediate value of m z = 10 10 M provides, as expected, re- 
sults that are intermediate between the two cases presented 
above: the observations are reproduced for r z = 500 kpc 
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Figure 2. The dependence of the relation between tHi and rCiv on the model parameters r z (the physical radius out to which galaxies 
enrich the IGM) and m z (the minimum halo mass that is capable of enriching the IGM). Ea ch pane l show s a different choice for m z 
and in each of the panels the yellow points with error bars show the observed relation of Aguir re et al.l (2005), with arrows representing 
upper limits. The curves show the simulation results, with the different colours representing the different choices for r z , and different 
line styles corresponding to different simulations. The metallicity, Z, required to match the normalisation of the observed data at the 
maximum value of tHi (for the maximum value of r z ) is shown at the top of each panel. The black curve in panel a is not visible because 
for ra z = 10 8 M© both r z = 250 kpc and r z = 500 kpc have fy = 1.00 and hence predict identical results. All coloured curves in each 
panel are scaled to the same metallicity. In order to match the observed median optical depth of C IV in regions with a low H I optical 
depth, it is necessary that low-mass haloes (ra z ^ 10 10-5 Mq) enrich the IGM out to distances > 10 2 kpc. 



but not for r z = 250 kpc and smaller. Considering now the 
extremely high-resolution simulation (L012N512; panel a), 
we find that if m z is as low as 10 8 M©, then the observed 
rc iv-thi relation can be matched even if r z is as small as 
62.5 kpc, but not for r z =31.25 kpc. 

Our results for the volume and mass filling factors of 
metals are summarised in Fig. [3l where we show the mass 
(upper panel) and volume (lower panel) filling fractions as a 
function of r z for different values of m z . Models that predict 
enough C IV at low th i are indicated with bold symbols and 
are connected with solid, thick lines. It is clear that in or- 
der to reproduce the observations, we require that the IGM 
is enriched primarily by low-mass (ratot < 1O 1O M0) haloes 
that are capable of driving gas out to distances > 10 2 kpc 
by z = 3. Allowing only haloes with masses ^ 10 10 - 5 M to 
enrich the IGM, results in metal distributions that fall far 
short of the observations for all sensible values of r z . Mod- 
els relying on (the progeni tors of) observed Lym an break 
galaxies (m to t ~ 10 12 M^ : lAdelberger et a"i1l2005h to do the 
enrichment, predict filling factors that are far too small to 
account for the observations and have not even been plot- 
ted. In our models, the volume filling factor of metals must 
be > 10% in order to reproduce the observations, indicat- 
ing that the observed rc iv - thi relation tells us robustly 
that a significant fraction of the volume and mass of the 
IGM is enriched by metals. On the other hand, a model vol- 
ume filling factor > 10% is not necessarily sufficient. Models 
that attain such large volume fractions by enriching the gas 
around high-mass haloes to very large distances often fail to 
reproduce the observations. We note that the values of the 
volume and mass filling fractions that we infer by comparing 
our models to the observations are approximate, even ignor- 
ing uncertainties associated with small-scale metal mixing. 
The filling fractions could change by factors of a few due to, 
for example, scatter in the metallicity around galaxies of a 
fixed mass and non-spherical bubbles. 



Finally, in Fig. 2] we show the effect of allowing haloes to 
enrich out a fixed multiple of their virial radii, r2oo- The yel- 
lo w data points with e rror bars show again the observation s 
of ISchave et al.1 (|2003T ) as presented in I Aguirre et al.1 (|2005h . 
and the curves show the results for the imposed metallic- 
ity distributions indicated in the legend. It is striking that 
even if all haloes that form stars are able to pollute the 
IGM as far out as four times their own virial radius, the low 
optical depth part of the rc iv - thi relation is not repro- 
duced. As dynamical processes (e.g. tidal or ram pressure 
stripping) are only expected to affect the metal distribution 
on scales < r2oo, this implies that some sort of ejective feed- 
back is necessary in order to pollute the low-density IGM. 
The conclusion that dynamical processes such as stripping 
ar e not, on their own, su fficient is consistent with the models 
of lAguirre et al. | (|2001ah as we l l as w ith the hydrodynamical 
simulations of IWiersma et al.1 (|201lh . 



4 DISCUSSION AND CONCLUSIONS 

By combining realistic cosmological density and tempera- 
ture distributions with toy models for the metal distribu- 
tion, in which only haloes more massive than a critical mass 
(m z ) are able to enrich a spherical region (of proper radius 
r z ) to a metallicity of Z, we have investigated which metal 
distributions can reproduce the observations of C IV asso- 
ciated with weak H I absorption (thi <C 10) as measured in 
quasar absorption spectra. 

The results presented in Sec. [3] imply that in order 
to match the observed median optical depth of C IV as a 
function of thi, we require that low-mass haloes, 771200 < 
10 10 M , are able to drive metals into the IGM, enriching 
the gas around them out to proper distances of > 10 2 kpc 
to metallicities of ^ 10 -3 Z©. 

We now verify that this scenario is physically possi- 
ble by comparing the required metal mass to an estimate 
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Figure 4. The 2 = 3 rCiv — tHi relation for models in which 
r z (the physical radius out to which galaxies enrich the IGM) 
is set to a fixed multiple of the halo virial radius. Yellow points 
with error bars represent the observations presented in Aguirre 
et al. (2005). Each curve represents a different imposed metal 
enrichment pattern, with different colours corresponding to dif- 
ferent values of ra z , and different line styles representing different 
choices for r z , as indicated in the legends. Even if all haloes can 
enrich the IGM out to four times their own virial radius, the ob- 
served rC iv — th i relation is not reproduced. This implies that 
dynamical effects alone are insufficient to explain the observed 
distribution of C IV. 



Figure 3. Fraction of the volume (/v; top panel) and gas mass 
(/m; bottom panel) of the simulation that is enriched with metals 
as a function of r z (the proper radius out to which galaxies enrich 
the IGM). Each colour represents a different choice for ra z (the 
minimum halo mass that is capable of enriching the IGM), and 
the diamonds (triangles) show results for L012N512 (L025N512). 
Thick curves and symbols show models that reproduce the rC iv — 
th i relation, thin, dotted lines show models that do not. We 
require fy > 0.1 and f m > 0.5 in order to match the observations. 



of the maximum allowed mass of carbon in the Universe. 
We can estimate the maximum allowed mass of carbon 
in the IGM from the stellar mass density at z = 3. Un- 
der the assumption that all of the stars formed in a sin- 
gle burst at t = 0, and assuming a Chabrier stellar ini- 
tial mass function and using the lifetimes, yields and su- 
pernova rates used in the simulation and summarised in 
IWiersma etHI ([2009b), the maximum allowed cosmic den- 
sity of carbon is 0.006p*, where p* is the cosmic density in 
stars and has been me asured to be p* = 10 7 ' 55 M /Mpc 3 
(Mar chesini et al.| [2009) at z — 3, giving us a maximum al- 
lowed total density in carbon of pc = 2.2 x 10 5 M /Mpc 3 . 
For the models that can successfully reproduce the observed 
optical depth distributions, the cosmic density of carbon 
varies from 8 x 10 3 M /Mpc 3 (ra z = 10 10 M ; r z = 250 kpc) 
to 2 x 10 4 M /Mpc 3 (ra z = 10 9 M ; r z = 500 kpc). In each 
case the total amount of carbon is an order of magnitude 



below the maximum allowed amount^, so the models are 
physically reasonable on these grounds. 

We thus conclude that in order to recover the observed 
median C IV optical depth in regions of low H I optical depth 
(thi rsj 1 — 10), we require that the galaxies in low-mass 
haloes (ra z < 1O 1O M ) enrich the IGM out to distances 

> 10 2 kpc. Galaxies residing in much higher mass haloes 
(ra z > 10 10 ' 5 M ) are too rare and too strongly clustered to 
contribute significantly to the enrichment of the low-density 
IGM. In every one of the models that is capable of repro- 
ducing the observations, the metal volume filling factor is 

> 10% and the gas mass fraction enriched with metals is 

> 50%. 

As discussed in detail in section 12.21 our models are 
not fully self-consistent in the sense that we have imposed 
metal distribution is post-processing. Although we found 
that post-processing hydro simulations without winds or 
with very strong winds led to identical conclusions, we can- 
not rule out the possitibity that future self- consistent models 
would yield different results. Unfortunately, current simula- 
tions suffer from large uncertainties due to their use of sub- 
grid recipes for the generation of winds and still lack the 
resolution to resolve outflows from low-mass galaxies in a 
representative volume. 

Two additional caveats must, however, be stressed. 
Firstly, the models presented here implicitly assume that 



2 Including a lognormal scatter in the metal optical depth at the 
level measured by Schav e~et al.l ([20031 ) would raise the metallicity 
by a factor of a few, which is still comfortably within the allowed 
limits. 
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the metals are we ll mixed on small sc ales. If, as suggested by 
observations (e.g. ISchave et al.l l2007). the intergalactic met- 
als are concentrated in metal-rich patches which together 
account for large covering factors, then the required filling 
factors could be smaller. In that case the filling factors de- 
rived here apply to the metal distribution smoothed over 
the scales that are somewhat smaller than the size of the 
bubbles, i.e. tens of kpc. 

Secondly, there exists considerable uncertainty in the 
spectral shape of the ionizing background, which leads to 
uncertainties t he fraction of carbon that exists as CIV. 
lAguirre et al.l (|2008l ) considered several models for the UV 
background, inclu ding some extreme one s, and found that 
only the fiducial Haardt &; Madaul (|2Q0lh model, which is 
the model used here, resulted in reasonable val ues for the rel- 
ative abundances inferred from observations. ISchave et al.l 
(2003) found that assuming a much harder (softer) spec- 
trum would flatten (steepen) the metallicity-density relation 
inferred from the observed tciv-thi relation, which implies 
that our models would predict steeper (flatter) tciv-thi re- 
lations than for our standard UV background, which would 
increase (decrease) the inferred filling factors. 



Our finding that the observations imply that the IGM 
was enriched by very low- mass galaxies is in agreement 
with a variety of the o retica l s tudies (e.g. lAguirre et aTl 



2001bl: Madau et all l200ll: 


ScannaDieco et al.l 120021: 


Thacker et al. 2002: Samui et al. 


120081: lOoDenheimer et al.l 


2009: Wiersma et al.1 hold). For example, bv modelling 



the propagat ion of galact ic winds in already completed 
simulations, lAguirre et all ([2001b) found that in order 
to explain the metallicities measured in the low column 
density part of the IGM, galaxies with baryonic masses 
> 10 8-5 M0 needed to launch winds with velocities of 
at least 200 — 300 km/s. The results pr esented here are 
also consistent with Wiersma et a D (|201Ch . who used fully 
self-consistent hydrodynamical simulations (including the 
simulations underlying our models) that massive haloes 
(> 10 12 M©) are unimportant for the enrichment of the 
diffuse IGM and that most of the metals that reside in the 
IGM at z = were ejected by a population of low-mass 
galaxies at high redshift. 

In contrast. IScannapieco et aTl (|2006l ) found that in or- 
der to match the strong clustering of C IV lines that they 
measured in their observations, the absorption needed to 
be generated primarily by gas that is strongly clustered 
around massive galaxies. Their best fit model required m z = 
10 12 M and r z = 500 kpc. Using these parameters in our 
model leads to such low volume (/v = 6 x 10 -3 ) and mass 
(/m = 0.03) filling factors that the corresponding curves 
would not appear on any of the plots presented here. We 
plan to explore the apparent tension between line clustering 
and optical depth statistics in future work. 

In summary, we have found that in order for a sim- 
ulated cosmological gas density field to reproduce the ob- 
served tciv-thi relation, we require that both the fractions 
of the volume and mass of the IGM that have been polluted 
with metals are substantial, with all our successful parame- 
ter choices yielding model volume filling factors greater than 
10% and model mass filling factors greater than 50%. The 
models favour metals being ejected from a population of 
low-mass (m to t < 10 10 M©) haloes at high redshift. 
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